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We investigate energy resolved electric current from various correlated host materials under out- 
of-equilibrium conditions. We find that, due to a combined effect of electron-electron interactions, 
non-equilibrium and multi-particle tunnelling, the energy resolved current is finite even above the 
Fermi edge of the host material. In most cases, the current density possesses a singularity at the 
Fermi level revealing novel manifestations of correlation effects in electron tunnelling. By means 
of the Keldysh non-equilibrium technique, the current density is calculated for one-dimensional 
interacting electron systems and for two-dimensional systems, both in the pure limit and in the 
presence of disorder. We then specialise to the field emission and provide a comprehensive theoretical 
study of this effect in carbon nanotubes. 
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I. INTRODUCTION 

At the early stages of the development of the quan- 
tum theory it became clear that electron tunnelling pro- 
cesses are of fundamental importance for condensed mat- 
ter physics [Q. Electron tunnelling became an essential 
concept in such fields as semiconductor physics, parti- 
cle transport in mesoscopic physics, field emission, and 
many others. In a multitude of setups, tunnelling pro- 
cesses have attracted attention over decades and continue 
to do so. 

For the purposes of this work, we shall visualise tun- 
nelling events as taking place between a metallic host 
material and a lead material (another metal, a semicon- 
ductor, or indeed vacuum as in the field emission effect). 
The host and the lead are separated from each other by 
a substantial potential barrier. The lead is supposed to 
contain a detector measuring the current density. The 
presence of a finite current immediately implies that one 
is dealing with a non-equilibrium (steady-state) system. 
Apart from the total current, the energy resolved current 
j{Lo), i. e. the amount of current in the energy window lu 
to Lo+du, is an important characteristics of the tunnelling 
process. The total current J is then the integral of j{uj) 
over all energies. In the field emission (FE) setup, the 
energy resolved current have been measured for different 
emitters over some 30 years (see 1^, and also references 
below) . We shall therefore adopt the terminology of the 
Gadzuk and Plummer's review article, (2|, and refer to 
as the total energy distribution or TED. We are not 
aware of any direct measurements of TEDs in tunnelling 
contacts. 

From the physical point of view, the most interesting 
aspect of the TED is the part of the spectrum above 
the Fermi edge, ui > Ep. At the leading order in the 
tunnelling amplitude, j{uj) is simply proportional to tun- 
nelling probability times the electron energy distribution 
function in the host, n{uj). The latter is identically zero 
(at zero temperature) above the Fermi edge due to the 
Luttinger theorem The high-energy tail is due to a 
combined effect of non-equilibrium, multi-particle tun- 



nelling processes, and mutual electron-electron interac- 
tions in the host material (for non-equilibrium multi- 
particle effects per se are not sufficient to smear the Fermi 
surface [Q). So, we shall also refer to this high-energy tail 
of the energy resolved current distribution as 'secondary' 
current, 'primary' current being the one below the Fermi 
edge (i. e. the one which is proportional to the tun- 
nelling probability). In the grand picture, the secondary 
current is akin to the Auger effect though it is, of course, 
a more complex phenomenon taking place not in an atom 
but in a fully interacting metallic host. Measurements of 
the secondary current (the trivial thermal broadening of 
the primary current being subtracted off) can thus pro- 
vide a valuable source of information about the electronic 
correlations in the host ||]. The effect was indeed 
first discovered experimentally - in FE measurements by 
Lea and Gomer Q. Theoretical analysis of this phe- 
nomenon was done by Gadzuk and Plummer soon after- 
wards . Following the pioneering paper by Fowler and 
Nordhcim these authors used the connection between 
the FE problem and the tunnelling problem and studied 
Boltzmann-like equations for the particle-hole balance in 
the low density approximation. To remove such restric- 
tions and to put the theory on the modern footing, so that 
it becomes applicable to strongly-correlated emitters, we 
embark in this paper on investigating the issue by em- 
ploying Keldysh diagram technique, appropriate for this 
non-equilibrium situation j|, |l^ . This method allows us 
to consistently write down series in the tunnelling am- 
plitude for all quantities of interest and for the TED in 
particular. 

Our motivation is in fact two-fold. On one hand, we 
think that the secondary current phenomenon has not 
received sufficient attention of theorists. The physics 
of the interplay between non-equilibrium multi-particle 
tunnelling and electron correlations is worth a deeper 
study. In particular, we investigate in this paper what 
aspects of the electron correlations are responsible for 
the current above the Fermi edge and in what setups, 
other than the FE, such current can occur. On the other 
hand, since the original work [0, Q there have been con- 
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siderable advances in the emitter technology. Perhaps 
the most important recent development is the usage of 
carbon nanotubes as field emitters. Carbon nanotubes 
display a remarkable array of electronic and mechani- 
cal properties and are of potential technological impor- 
tance . Electron transport in these systems has been 
thoroughly investigated. While single-wall nanotubes 
(SWNTs) exhibit one-dimensional (ID) Luttinger liquid 
(LL) type transport prop erties and are pretty well un- 
derstood, the theory J14 |l3| , and experiment being 
in agreement multi-wall nanotubes (MWNTs) 

are more complex systems subject to intensive current 
debate. MWNTs are composed from many (at least ten) 
concentric graphite shells. Current experiments on them 
are consistent with two-dimensional (2D) diffusion with 
characteristic weak localisation features and zero- 
bias anomalies [Esl EoL EOl 21 . What concerns us here 



is that, apart from other uses, carbon nanotubes are ex- 
pected to act as field emitters in high-resolution displays 
and cathode tubes ^ . While there have been several 
experimental investigations of the FE from carbon nan- 
otubes (see |2^ and the main text for more references), 
both from SWNTs and MWNTs, the relevant theory has 
been lacking. So, the second leg of our motivation is 
to discuss existing FE experiments on carbon nanotubes 
and make further theoretical predictions for these sys- 
tems. 

Having in mind applications to carbon nanotubes and 
taking into account the fact that FE usually occurs from 
a tip of the tube (both for SWNTs and MWNTs) we 
narrow the following considerations from a number of 
imaginable setups to an appropriate tip geometry. Apart 
from this restriction, wc intent to advance in this paper a 
general discussion of non-equilibrium multi-particle tun- 
nelling from strongly-correlated ID and 2D hosts paying 
special attention to carbon nanotubes. Some of the re- 
sults on ID emitters (applicable to SWNTs) have been 
announced in our recent letter 

The paper is organised as follows. In the next Sec- 
tion we present some qualitative considerations concern- 
ing the physical nature of the TED. Section III contains 



general (model-independent) results. We identify the rel- 
evant Keldysh diagrams contributing to the TED above 
the Fermi edge and then perform a spectral analysis of 
the involved correlation functions. A simple application 
of the developed theory is contained in Section IV, where 



we treat an electron system with local interactions con- 
fined to the vicinity of the tunnelling point. In Section 
^ we analyse the TED of particles tunnelling from LLs. 
In subsequent Sections Vl and VII, the same problem is 
studied for correlated 2D electron systems, respectively 
in the pure limit and in the presence of a disorder po- 
tential. In Section VIII we specialise to the case of field 
emission from carbon nanotubes. Finally, in Section IX 
we discuss a more sophisticated two-stage tunnelling via 
a localised state. Summary and conclusions Section com- 
pletes the paper. 



II. PHYSICAL PICTURE 

Before proceeding with calculations, let us elaborate 
on qualitative origins of the TED high-energy tails. The 
simplest setting to start with involves two noninteract- 
ing electrodes with a tunnelling contact between them, 
see Fig. |l[ Applying a finite voltage leads to a nonzero 
current through the junction. Neglecting the charging 
effects and (for now) under the assumption of constant 
electron densities of states and transmission coefficient, 
the current is proportional to the applied voltage V. On 
the right side of the contact, where the chemical potential 
is supposed to be lower than on the left one, the current 
is carried by the tunnelled electrons. On the left side of 
the contact the current is carried by the holes moving in 
the opposite direction away from the contact. As long as 
the system is noninteracting, the TED of the electrons 
that tunnelled out is uniform in the window between Ep 
and Ep — V and is zero outside (this is evident but can be 
confirmed by simple calculation in the spirit of jj] , which 
we omit). From now on we set e = 1. Ep is the Fermi 
level of the left electrode. 

This picture changes drastically if interactions between 
the electrons are switched on. We restrict our consider- 
ations to the case when only the left contact - the host 
- is correlated (see also discussion in the next Section). 
Then the difference between the actual energy distribu- 
tion function in the lead and the noninteracting distribu- 
tion function gives the TED of particles tunnelling out 
of the host. The holes left behind by electrons tunnelling 
out below the Fermi energy ('primary' electrons) expe- 
rience scattering from the electron sea, thereby creating 
electron-hole pairs. Contrary to the problem of hot elec- 
tron relaxation (see e. g. p5| , and references therein), we 
are dealing here with a flow of hot holes in a steady state. 
Also, we are interested in a different object: the energy 
distribution functions rather than momentum distribu- 
tion functions (the latter do, of course, have a non-zero 
tail above the Fermi momentum solely due to correla- 
tions, without a need to include tunnelling [^). The 
emerging 'secondary' electrons can also be carried over 
to the lead by a successive tunnelling process, as shown 
in Fig. 1^. This phenomenon can be regarded as con- 
densed matter analogy of the Auger process known from 
the atomic physics. Since the energy of the holes, mea- 
sured from the Fermi edge Ep, can not exceed V, the 
upper limit for the energy the secondary electrons can 
acquire in the electron-hole pair creation process is given 
by V as well. At higher orders in the tunnelling prob- 
ability this threshold is increased so that the secondary 
electrons can, in fact, become more energetic than V. 
In this case they ought to be successively scattered ei- 
ther from another hot hole or from another hot electron. 
This is, however, a process of higher order not only in 
the tunnelling but also in the interaction. Therefore sec- 
ondary electrons with energies between {n—l)V and nV 
emerge in the processes of 27T.th order with respect to the 
tunnelling amplitude as well as the interaction constant. 
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FIG. 1: Tunnelling junction biased by finite voltage V. The 
left electrode is the interacting one while the right one is un- 
correlated. Tunnelling electrons with energies higher than Ef 
are due to hot hole (denoted by Pi) scattering with creation 
of electron- hole pairs e — pc (see inset). 



Obviously, in equilibrium V = all the high energy tails 
vanish, thereby restoring the Fermi edge in accordance 
with the Luttinger theorem. So, non-trivial interaction 
effects are encoded in the high-energy behaviour of the 
non-equilibrium TED. 

Theoretically, the most interesting limiting case is the 
behaviour of the TED just above the Fermi edge, when 
certain universality can be expected. Out of the consider- 
ations of Ref. a divergent TED emerged, with (simpli- 
fying matters) the singularity approximately of the form 
j(aj) ~ — Ep)- This turned out to be roughly con- 
sistent with measurements at the time As detailed 
in the following Sections, the present study supports the 
view that the limiting form of the TED strongly depends 
on the nature of the host material and the geometry of 
the setup (within a given material type though, some uni- 
versality does set in, so for point-contacted LLs we find 
a power-law behaviour, etc.). 



III. STATEMENT OF THE PROBLEM AND 
GENERAL RESULTS 

We now formalise the problem by writing down the 
relevant tunnelling Hamiltonian: 

H = H[c] + H[i>] + J [ij^O)c{0) + c^O)iP{0)] . (1) 

Here 7 is the tunnelling amplitude and c and ijj are the 
annihilation operators for the electrons in the lead and 
in the host, respectively. The unperturbed part of the 
Hamiltonian Hq = H[c] + H[tp] describes two decoupled 
electron systems at different chemical potentials /x^ — 
fic^V>0. 

A clarification on the following points is in order. 



(i) As we only want to consider a tip geometry, we have 
explicitly assumed that the tunnelling occurs only locally 
at the location of the tip: x = 0. In reality, of course, 
there is a small area over which the tunnelling takes place 
(strictly speaking, a real-space integral is required in the 
tunnelling term in Eq.(|l|)). As this has no qualitative 
influence, we shall keep writing simple as long as we can. 
The locality assumption is natural for carbon nanotubes 
because of the very shape of these objects. It is however 
also justified for most bulk interfaces and field emitters 
where, because of the roughness of the surface and the 
pronounced exponential dependence of the penetration 
coefficient on the distance between the electrodes, the 
main contribution to the current comes from only few 
points between the electrodes. 

(ii) Related to the above is the question of the energy 
(momentum) dependence of the tunnelling amplitude. In 
real systems it is energy dependent. So, for the EE setup 
the relevant energy scale is determined by a combination 
of the work function and the applied field. It is impor- 
tant to keep in mind, however, that as the tunnelling 
amplitude is a single electron property, regarding its en- 
ergy dependence there is no special significance to the 
Fermi edge. Therefore, when addressing observables de- 
termined by scattering processes taking place close to the 
Fermi surface, it is quite safe to neglect the energy de- 
pendence of the tunnelling amplitude and replace it by 
a constant, 7 (7^ being proportional to the transmission 
coefficient of the potential barrier at the Fermi energy: 
'D{Ep).) In the literature on the subject the tunnelling 
term is often quoted in the momentum representation 
in the form jpkipj^Ck plus conjugate with some unspec- 
ified matrix elements 7pfc. We take the view that, for 
calculating quantities related to the Fermi surface, this 
would only complicate formulae without conceptual gain. 
There will be, however, instances in the following when 
the energy dependence of the tunnelling amplitude is im- 
portant, like when there is a localised state or when we 
discuss generalisations of the Fowler-Nordheim relations 
for nanotubes. In these cases (Sections VIH and IX) we 
shall take the relevant energy dependence fully into ac- 
count. 

(iii) Throughout this paper we assume that the lead 
is non-interacting. Indeed one could not otherwise sepa- 
rate the effect of correlations occurring in the host from 
those occurring in the lead, which would seriously ham- 
per meaningful interpretation of measurements on such 
systems. This assumption is justified for the FE setup, 
save for the Boersh effect, which is not believed to be 
important for carbon nanotubes (see and ||]). For 
a general tunnelling junction setup it would not be cor- 
rect to a priori neglect correlations in the lead. However, 
with care it is possible to realise a reasonable setup in- 
volving, for instance, a nanotube in contact with a low 
carrier density semiconductor or a quantum wire opening 
up into a higher dimensional lead preferably screened by 
a nearby gate (in fact most conductance measurements 
on quantum wires are nowadays interpreted in terms of 
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FIG. 2: Schematic representation of the Keldysh contour. 
The lower branch is time-ordered, while the upper branch is 
anti-time-ordered. 



a junction with a non- interacting lead p9[). 

In the bulk of this paper we shall be investigating 
the secondary current j{uj > 0) (we set Ep ~ from 
now on). It is a plausible statement that the secondary 
current should be proportional to the high-energy tail 
nlco > 0) of the electron energy distribution function 
inside the emitter tip (calculated to all orders in the 
tunnelling amplitude). Indeed the latter quantity corre- 
sponds to the amount of electrons available for tunnelling 
at a given energy lo, while the proportionality factor con- 
tains (unessential) information about how these electrons 
are then carried over to the detector. In the first part of 
this Section we establish this statement. 

Since we are dealing with non-equilibrium phenom- 
ena we resort to the Keldysh formalism. We denote by 
G{x,x';t — t') the generalised Green's function for elec- 
trons in the lead, 



G{x,x']t - t') = -i{Tc[c{x,t)c'{x' ,t') Sc])a 



(2) 



the Green's function in the host, g, being defined by an 
analogous formula. Here Tc stands for the ordering op- 
eration along the Keldysh contour, shown in Fig. ^, and 
X corresponds to a set of coordinates specifying the elec- 
tron states in the lead. (It should be understood as a 
distance from the emitter tip plus possibly a transverse 
channel index which we suppress as it plays no impor- 
tant part in the following). The average in Eq. (^) is 
taken over the ground-states of the unperturbed Hamil- 
tonian Hq and the contour ^c-matrix is responsible for 
tunnelling events: 



Sc 



c-exp 



Z7 / dt['ilj\t)c{t) + c\t)^{t)] 



Note that because the system is not translationally in- 
variant, all Green's functions depend on both coordi- 
nates. However, all functions still only depend on the 
time differences because the system is in steady state. 
The time integration in Eq. (g) is along the contour C . 
Keldysh disentanglement of the time variables results in 
expressions with integrations only along the real axis. 
Then, four different Green's functions emerge in accor- 
dance with four possibilities to arrange the times t and t' 
along the contour. This placement of the time variables 
is refiected in additional superscripts of the Green's func- 
tions. For instance, if the time t lies on the time-ordered 



FIG. 3: Schematic representation of the equation (^. Solid 
line denotes the electron's Green's function of the lead and 
the dashed ones stand for particles in the host. Thick lines 
they correspond to exact Green's functions to all orders in 
tunnelling and interactions while thin lines represent Green's 
functions with tunnelling neglected. 



part of the contour and t' on the anti-time-ordered one, 
see Fig. |^, evaluation of Eq. (||) at the zeroth order in 
the tunnelling (indicated by a subscript 0) yields: 

G^+{x,x'-t-t') = -i{Tc[c{x,t)c^{x',t')])o 
= i{c\x' ,t')c{x,t))Q . 

Its counterpart with an interchanged orientation of time 
variables {t on the C+ and t' on the C_) is 

G'(j'~(x, x' ■. t — t') — —i{c{x, t)c\x' , t'))o . 

Green's function in which both time variables lie on the 
same side of the contour are the usual time-ordered and 
anti-time-ordered ones: 

Go"(a;,a;';i-i') = ~i{T[c{x,t)c^ {x ,t')])o , 
G++{x,x';t^t') = -i{f[c{x,t)cHx',t')])o, 

where T denotes the anti-time-ordering operation. 

The local electron energy distribution function N{x,ll!) 
in the lead, which we shall also call TED by abuse of 
terminology, is given by the defining relation with one of 
the Keldysh Green's functions (sec c. g. |10|): 



N{x,uj) — ~iG '^{x,x;uj) 



(3) 



The lead being non-interacting, the TED of the tun- 
nelling particles is, up to a pre-factor, given by above 
Green's function after subtracting off the equilibrium dis- 
tribution function. Indeed, by examining the perturba- 
tive expansion in the tunnelling amplitude 7 and keeping 
in mind the fact that there is no correlations in the lead, 
one can easily establish the following important identity. 



G{x,x';t-t') = G(i{x,x';t-t') 
+ 7 



(4) 



2 ' ' dt"dt"'Go{x,Q-t-t")9{'^A,t" -t'") 
c c 



X Go(0,a:';t"'-t') . 

The corresponding diagram is shown in Fig. |^. The func- 
tion (7(0, Q;t — t') appearing in this relation is the exact 
one, so the formula is valid for arbitrary interactions in 
the host and to all orders in the tunnelling amplitude. 

By disentangling the Keldysh indices and changing 
over to the energy representation, we extract the Green's 
function of interest from the general expression, (^, 

G-+{x, x'; uj) = Go +(a;, x'; ou) (5) 
+ 7' ^ (v;j)Go-''(x,0;c.)g^^(0,0;c.)G^o+(0,x';c.). 
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Plugging this expression into Eq. (||) one obtains the com- 
plete TED at all energies. Let us take a closer look at this 
relation. All terms on the right-hand-side, apart of the 
one with = (—,4-), contain the imperturbed Gq^ 

function. The latter is, in turn, proportional to the TED 
in equilibrium. Therefore, no matter what the Green's 
functions of the interacting fermions look like, all these 
contributions represent the TED below the Fermi energy 
(where the first term is the dominant contribution). The 
term in Eq. (|^) with (i,j) = (— ,+) is not constrained 
in such way and can, in principle, contribute above the 
Fermi edge. Therefore the high-energy part of the TED 
is given by 

N>{x;lu) = -ie{uj)G-+ix,x;uj) (6) 
= (x,0;c^)g;+(0,0;c^)G++(0,x;c^) , 

where Q{uj) is the Heaviside step function and 
g^'^{0, 0; uj) stands for the high energy part of the corre- 
sponding Green's function. 

Similar properties can be established for the energy 
resolved current j{uj), when we have 

1 f 

j{uj) = — / dkvkN(k;Lj) 

TT J 

where Vk is the velocity of the particle with wave number 
k. Making use of Eqs. (||) and (^) for the high-energy 
part of the current one obtains, 

.7>H = -z^.9> + (0,0;o.) (7) 

TT 

X J dkkj d(a; - a;') e*'=(^-^') 

where the product in brackets is actually translationally 
invariant (only depends on a: — x'), reflecting the fact that 
the excess particles only travel in one direction, away 
from the contact, in the (non-interacting) lead. 

We now observe that both Eq. ^ and Eq. (^ are pro- 
portional to the high-energy part of the TED of particles 
at the tip of the interacting host, 

n>(c^) = -^5-+(0,0;a;), (8) 

thereby proving the statement put forward at the be- 
ginning of this Section. We stress again that ny{uj) is 
supposed to be exact both with respect to the interac- 
tion and the tunnelling. Indeed, as one can easily see, 
this quantity is zero for a noninteracting system, even if 
tunnelling is taken into account to all orders. It is still 
zero for interacting hosts if tunnelling is neglected or if 
the system is in equilibrium (Luttinger theorem). Situ- 
ation changes dramatically in the case of an interacting 
host, finite tunnelling and finite bias voltage. In the rest 
of this Section we shall explore what statements can be 
made about the TED without assuming a concrete model 
for the host material (other than that there is a two-body 
interaction present). 



a) b) 




+ + 



FIG. 4: a) The only second order diagram contributing to 
the TED above the Fermi energy. Solid lines correspond to 
host electrons and the dashed one to those in the lead while 
crosses denote the tunnelling vertices. Wiggly lines represent 
interactions, b) This diagram has to be inserted into a) at 
points denoted by circles in order to obtain high-order dia- 
gram for the TED. 

As mentioned in the Introduction, there is a scattering 
process that allows for creation of electrons above the 
Fermi energy. Let us start with perturbative analysis. 
At the second order in the interaction, the correspond- 
ing diagram for the Green's function can be constructed 
by annihilating the particles in the same way they were 
created. This essentially converts the scattering ampli- 
tude for a given process into the corresponding proba- 
bility. The result is shown in Fig. ^ a). Now we need 
to establish the way the vertices can be decorated with 
the Keldysh indices. Trivially the outmost points should 
have index — on the outward leg and + on the inward 
leg in accordance with the type of function ( — h) we are 
calculating. Furthermore, there is only one possibility to 
assign indices to interaction vertices. In order to obtain a 
contribution above the Fermi energy the inward Green's 
function should be of or anti-time ordered type while 

the outward one has to carry indices, because all 

other possibilities contain at least one Gq^ factor which 
forces the diagram to vanish above Ep. We are left with 
two cross vertices of the inserted Green's function. Since 
we do not have a priori knowledge about them, we have 
checked all possibilities explicitly. It turns out that the 
only possibility is to insert the function G^^ because in 
all other cases the diagram vanishes for ui > 0. Hence 
the only second-order diagram contributing to the TED 
above the Fermi energy is the one shown in Fig. |^ a). 

Let us now consider all orders in the interaction but 
remain at the second order in the tunnelling amplitude 
thereby also formally justifying the above choice of the 
Keldysh indices. The way to proceed is to analyse the 
relevant Green's functions in the Lehmann-type spec- 
tral representation. In the time domain the local (tip) 
Green's function of interest g ''(ti, is given by 

g-+{h,t2) = i{i>Ht2mti)) = -^{TcmtT)^H4)]) , 

where the superscripts indicate that the time variables ti 
and t2 lie on the time-ordered (T) and anti-time-ordered 
(T) parts of the Keldysh contour, respectively. Perform- 
ing the straightforward S'-matrix expansion in powers of 
the tunnelling amplitude 7 one obtains at the lowest non- 
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vanishing order (second order in 7): 



ij = ± 



(9) 



where /Cj,~'"(i2, ^i; ^3, ^4) are the foUowing four-point cor- 
relation functions: 



lC-_X{tut2-MM) = {f[i^{U)i^\h)]T[^{h)^\h)])o, 
1Czt{tiM;hM) = {^\t2m{ti)^\h)ij{U)])a, 
lC-X{hM;hM) = {T[i^Kt2W{hW{ti)]^{h))Qm 

We shaU first find the spectral representation of the 
term containing /C^^^, 

X / dt3di4e-"(*^'**)/C;t(ii,t2;i3,t4) 

(11) 

Green's function in the 



-^72 / dep,(e)p-+(ii,t2;e), 

-V 



where we have used the bare 
lead, 



G+-{e)^^i2TiQ{e + V)p,[e), 



(12) 



containing the local density of states Pc{^) in the lead. 
Partial Fourier transform of the correlation function ap- 
pearing in the above formula is defined by 

P+ + (ii,t2;e) = j dhdUe-'<'^-'^^lC+t{ti,t2-MM)- 

According to definition (^0|), this correlation function 
contains two different time ordcrings. Writing them down 
explicitly and inserting between every two ip operators a 
complete set of exact states, it is possible to perform all 
time integrations. This procedure is a straightforward 
generalisation of the standard one for the equilibrium 
case The result is 

P^-+(ii,t2;e)-^e-(^-+^)(*-*^)|SW(6)p, (13) 



with B\j. (e) defined by 



Ef,- E^ + e + iO E,y + e-iQ 



where Greek indices count all possible excited states of 
the system with energies £^!y,A,/i and a^j/ stand for matrix 
elements of the operator ^ , a^jy = {fi\'ijj\h') . In order 
to obtain the actual correction to the TED we plug the 
two last equations back into Eq. (O) and compute the 



Fourier transform of the latter with respect to the time 
difference ti — t2, 



(5n(l)(w) = -^ J d{h-t2)e'^^*'-''U^^^g-+{ti,t2) 



Obviously, all iJ^'s are larger than the ground state en- 
ergy Eq (which we have set to zero). Therefore the upper 
boundary for uj is given by V. 

The remaining three terms in Eq. (^ can be treated 
in a similar manner. Repeating steps leading to (p^), 
we obtain for the term containing the second four-point 
correlation function in Eq. (10): 



PI+{h,t2;e)^J2^ 



with 



Ef,-e + iO E^- E^-e- M 



the corresponding contribution to the TED being 

5(2)n(c.) = 2717^ - - ^)\B'^^\V)? ■ (14) 



Contrary to Eq. (|lj), this term is bounded by — from 
above and hence does not contribute to the high-energy 
part of the TED. One can use an even simpler argument 
in order to show that the last two four-point correlation 
functions in ( p^ do not contribute above the Fermi edge 
either. By inserting only one complete set of states be- 
tween the time-ordered operators (i. e. at the break of the 
time-ordering) one obtains (for simplicity we set ti = 0) 



Pr+(0,t2;e) 



dt^ dt^ e 



-ie(t3-t4) 



E aS.e-^^*MHT[V^(0)^nt3)V^(t4)] |0) 



for the third term. The corresponding correction to the 
TED is then given by 



J(3)n(c^) =Z72EaS.^(-^- -^) [deGo 

dHdtMT\^{^W{^^)^{u)m ■ 



Since all EJ'^ are always positive this expression is 
nonzero only for negative energies w < 0. The same 
is true for the last four-point correlation function /C^J. 

The above approach is quite general (i. e. valid for all 
kinds of interacting host materials) but it is also incon- 
clusive in the sense that the actual energy dependence 
of the TED is determined by spectral weights encoded 
in the structure of matrix elements a^i^ that is different 
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FIG. 5: Generic diagram giving contributing to the TED 
above the Fermi energy. The four-point correlation function 
fCT'!i{ti,t2;t3,t4) is represented by the shaded square. 



for different hosts. One positive result of the spectral 
method, however, is that we have identified the Keldysh 
four-point correlation function responsible for the sec- 
ondary current effect: IC^_, see Fig. |. For some models 
we shall calculate this correlation function exactly to all 
orders in the interaction, for other models we'll resort to 
perturbative expansions: note that diagram a) of Fig. ^ 
is a special case (second-oder expansion) of the general 
Fig. H term. 

Let us conclude the present Section by making one 
more observation of a general character. Returning back 
to the real time representation of the lead Green's func- 
tion ( p^ and doing simple manipulations with the time 
integrations, we arrive at the alternative representation 
for the TED: 

CxD OO OC 



t + ia 



— oo — oo 



X {T[^j{Ti)^{0)m{t + Ti)^{t + n + r2)])o , (15) 

where a stands for the short-time (high-energy) cutoff 
inversely proportional to the conductance band- width D. 
For the sake of simplicity, we have assumed a constant 
density of states in the lead but this assumption is not 
crucial. We now observe that if the T-ordering operation 



in the above formula were dropped, the expression would 
vanish. Indeed, inserting some complete sets one finds 
that then the integrand in the above, proportional to 



{Ol^Pit + Tl)\v){v\ij{t + n + T2)\0) 



(16) 



becomes an analytic function of the time variable T2 in 
the upper half plane and hence has vanishing Fourier 
transform for lo > 0. By re- arranging time integrations 
in a different way, one can easily show that the same 
statement is true about the T-ordering operation. Sub- 
tracting off the unordered correlation function we arrive 
at the following remarkable representation for the TED: 



oo oo 



-oo 

xn{Ti,0;t + n,t + Ti+T2) , (17) 

with 

n{hM-MM) = {{i^\h),ij\t2)}{^lj{h),ij{U)}) , (18) 

where {., .} stands for the anti-commutator. The advan- 
tage of this representation is that it is explicitly vanishing 
in the non-interacting case as the field operators are then 
anti-commuting at all times. The latter does not take 
place in interacting systems. Generally, anti-commuting 
interacting Fermi operators at different times would re- 
sult in a complicated object. However, there are systems 
for which the so-called braiding relations hold: 



V'(ii)V'(i2) = e''"(*i-*^VXi2)^(ii). 



(19) 



The exact shape of the phase function ip{t) depends on 
the system in question but it is usually proportional to 
the interaction constant. The braiding relation reflects 
the time evolution of an anti-commutators in particular 
models and can be used to simplify dl^) to 



J 



7^(Tl, 0; t + n, i + n + T2) = (1 + e-*'^("i')(l + e-'''^(-"^')(V'^(Ti)V't(0)V'(i + Ti)i>{t + n + T2)) 



(20) 



from which formula it is clear that the expansion of n> (w) 
in the interaction constant starts (at least) at the second 
order. In the next Section we present two models where 
such braiding relations exist. 



IV. MODELS WITH LOCAL INTERACTIONS 



We start with the simplest (toy) model we could think 
of: the interactions are only present at one point - the tip 
of the emitter. The host Hamiltonian entering Eq. (||) is 
then of the form 



In the following four Sections we apply the formalism 
developed here to various physical systems. 



(21) 
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with the Hubbard on-site repulsion for the interaction 
term 



i?, = C/o^|(0)Vl(0)Vi(0)VT(0) 



(22) 



We take a half-infinite electron system with linear dis- 
persion relation for the left electrode 



dx 



and the same for the lead (but, of course, at chemical 
potential of ~V). Here V's_r(l) stands for the annihilation 
operator of the right (left)-moving electron species with 
spin s. To simplify the formulae we have set hvp — 1- 

At first sight such model seems unphysical. However, 
it can be regarded as a special case of a correlated quan- 
tum dot contacted by two noninteracting electrodes, see 
e. g. Ref. where the contact to one of them is nearly 
perfect while the contact to the other one is weak. One 
should then decorate the electron field operator with a 
tunnelling channel (transverse quantisation) index. Its 
presence does not affect qualitative results, so we drop 
the channel index in this Section. 

For the half- infinite host, the bare Green's functions 
have image structure of the form 

9so{x, x'; uj) = go{x - x'; uj) ~ go{x + x'; uj) , (23) 

and are diagonal and s — > — s symmetric in the spin 
space. The Green's functions appearing on the right- 
hand-side in the above formula are the translationally 
invariant (bulk) ones for the corresponding spinless in- 
finite system. The simplest Green's functions are the 
retarded and advanced ones, 

g^ix; uo) = [g^{x- lu)]* = -»e^('=-+")l-l (24) 

to which the diagonal Keldysh functions are related via 

g^' {x-uj) = Q{uj)g^{x-^) + Q{-u)g^{x-uo) (25) 

and gQ^{x]Uj) — —[gQ^{x]Uj)\*, while the off-diagonal 
Keldysh functions are: 



+-(-+) 



(x; w) = i20(±a;) cos {kp-\'Uj)x 



Plugging these expressions into Eq. (^^ one obtains all 
Keldysh components for the half-infinite system. How- 
ever, using them directly to evaluate g7>^(0, 0; to) for elec- 
trons with spin orientation s causes a technical problem 
as at the origin all wave functions of the half-infinite sys- 
tem vanish and so do the correlation functions. The rea- 
son for that is the open boundary condition imposed on 
the wave functions. Therefore we are in need of a reg- 
ularisation. We now assume that the underlying lattice 
model has a lattice constant oq so that the tunnelling 
occurs between points with the spatial coordinate x < 



in the left subsystem as well as in the right one. Obvi- 
ously, such regularisation does not influence the physics 
but only creates some non-universal numerical factors in 
the TED. In this regularisation all zeros in Eq. (^ are 
understood to be substituted by oq. Suppressing this 
argument in the formulas, we obtain the following ex- 
pression for the second-order correction to the local (tip) 
Green's function (diagram Fig. ^ with an appropriate ar- 
rangement of the spin indices): 



--lHulQ{u;)g++{u,)g^-{co) 



X / deide2g„ +(w - ei) 



X 5o"+(ei+e2)go"-(e2)go++(e2)G+-(e2) (26) 

where the local Green's functions in the lead (G"s) are 
understood to be the same as in the host but with energy 
measured from — V. To evaluate this diagram, it is useful 
to remember that 



-4 sin^ [{kp + ijj)ao] 



As we restrict our considerations to energy scales much 
smaller than the Fermi energy, we can regard these prod- 
ucts essentially constant. The same is justified for other 
Green's functions except, of course, the step functions 
appearing in front of the off-diagonal Keldysh functions. 
Then, calculating energy integrals in Eq. ([2^), we obtain 
the following result for the TED: 

n,y{cu) = -^57>+(a;) = Co 7' &{iu)Q{V - - Vf/2 

(with the non-universal constant Co — 4^ sin^°(fci?ao)). 

Note that in accordance with the discussion in Sec- 
tion II no particles can have energy exceeding V at the 
second order in the interaction constant. Such particles, 
however, appear if one takes into account processes of 
higher order. The corresponding diagrams can easily be 
constructed, see Fig. ^. So, at the fourth order in the 
Hubbard U we found 

g^+{uj) ^iC^j'^ C/q e(tj)e(2F - uj){uj - 2T/)V24 

(Cq is again a non- universal numerical pre-factor). It is 
not difficult to calculate the particle spectrum for arbi- 
trary iv. In the energy window between V{n—1) and Vn 
{n = 1,2,...) it is dominated by the term of the order 
[7f/o]^" and shows {nV — loY^ decay. 

Another interesting local model is the local phonon 
model, where the electron density operator is coupled to a 
local oscillator at the tip a; = 0. In the chiral formulation 
the Hamiltonian of the problem is 



H = i J dxilj^ {x)dxij-'{x) 

+ av/(o)^a(o)(&^ + &) 



(27) 



where b^ and b are the creation and annihilation opera- 
tors of the local phonon. The chirality of fermions takes 
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care of the reflecting boundary condition at the origin 
(see e. g. Ref. where the tunneUing processes take 
place. In the anti-adiabatic approximation (justified for 
high-frequency phonons, Q 3> A/ao), this Hamihonian 
can be solved via the Frohlich transformation 

H' = , (28) 

with 

^^t(o)^(0)(6t_5). (29) 



Simple transformations of the field operators, which we 
omit here, show that braiding relations are satisfied for 
this particular model and the braiding phase (p defined 
in Eq. ( ^ is equal to: 

<^(t) = 7r+i (^A.^ sinm. (30) 

Using this result when evaluating Eqs. (|l5|) and (|l^) we 
obtain for the TED: 



oo oc 



n>(Lj) ^ ~2tJ^ / dty ^ / dTi / dT■2e»^("l+^^)(l-e*(V^ao)^si„[nrl]/2^.^_g-^(A/nao)=sin[^^r.]/2^ 

[t + iaY 





T1T2 



(ti + T2 + t + ia){Ti + < + ia){T2 + t + ia) 

V 



Expanding in powers of A/f2ao the evaluation of these 
integrals yields the main contribution to the TED of the 
form 



n>(w) 



As expected, the emerging spectrum has a resonant char- 
acter having a sharp pole at the oscillator frequency fJ. 
Under appropriate conditions, the Frohlich transforma- 
tion can also be used for solving the bulk electron-phonon 
interaction where the braiding relations still persist. We 
shall not pursue this issue further in this paper. 



V. LUTTINGER LIQUID MODEL 

Staying with ID hosts, we now move on to the next 
level of difficulty and discuss the bulk interactions. The 



simplest model here is the spinless Luttinger liquid model 
^3| . Subject to minor modifications (see below) the 
LL results will also be applicable to quantum wires and 
SWNTs. The interacting term in the host Hamiltonian 
then is 



H, = / dxdyxp\x)i;\y)U{x-y)ij{y)i;{x). (32) 



where U {x) is the interaction potential. 

Although this model can be exactly solved, it is in- 
structive to start with the perturbative expansion. Dia- 
gram a) in Fig. ^ still represents the only non-vanishing 
contribution to the TED above the Fermi edge. Com- 
pared to Eq. (^^, the analytic expression for this dia- 
gram is slightly more complicated since all participating 
Green's functions acquire spatial dependence. Assuming 
local interaction, U{x — y) UoS{x — y), we write 



g> + (w) = UqQ{lli) / dxidx2go{x2,ao;uj)g„ {ao,xi;uj) j deidt2 9^ {xi,X2]Uj - ni) 

X 9o^{xi,X2;ei +e2)go^{xi,aQ;e2)g^^{ao,X2;e2)G^^{ao,aQ;e2) . (33) 

We use the same regularisation as in Section ^ (that is why the lower boundary of the z-integration is the lattice 
constant ag). Substituting the expressions for the corresponding bare Green's functions into this formula one obtains 
the following result: 

g^+{cj) ^ i-f^U^SCoQiioMV^u;) f dci f ' dea |if(ei, 63, c^)P , (34) 

Jul J-v 

where Co is a non-universal numerical constant and the function H is defined as 

H{ei,e2,io)^l- f dxe^('"-^^)"cos[(2ei + e2-u;)x] ^ " 



2 Jq Ai ei{ei + 62 - lo) 



10 



We are ignoring terms containing rapidly oscillating inte- 
grands which are negligible on energy scales much smaller 
than Ep. 

Performing the last energy integrals we obtain for the 
TED above the Fermi edge 

n>H = eHe(v^ - c.)Co{-^ + Fo(^) } , 

where Fo{x) is a function regular at x = and vanishing 
for X > 1. While the TED vanishes smoothly in the 
vicinity of w = there is a sharp singularity towards 
the Fermi edge. 

Since the Luttinger liquid model is solvable for arbi- 
trary interaction strength, |33j, it is possible to go be- 
yond the perturbative expansion. In particular, the four- 
point correlation function in Eq. ( |l5| ) can be calculated 
exactly. We remind the reader that the system under 
consideration is a half-infinite spinless LL with an open 
boundary across which tunnelling processes occur. Us- 
ing the standard bosonization scheme for open boundary 
LLs (technical details can be found in the literature, see 
|l|, ||), we write 

^{x = 0,t) = (27rao)-i/2 exp[i(j){x = 0, t)/^ (35) 

where g is the LL parameter, defined by g = (1 + 



4!7o/7r) The Gaussian chiral Bose field 0(x,t) has 

been rescaled and is governed by the LL Hamiltonian 

Hll[4>] = dx{d^4>f. (36) 

The Bose field is periodic in 2L, where L is the system 
size. In order to obtain the correct analytic properties 
of correlation functions we first work with a system of 
a finite length L, which makes the energy quantisation 
equal to eo = tt/X, and then send L to infinity. Using 
representation (|3^) we obtain the four-point correlation 
function in question: 

lClt{tiM;t3.,U) = (27r6)-2sgn(t2 - h)sgn{ti - U) 
r F{\h-U\)F{\h-U\)FHO) 
VF{t2-h)F{t2~U)F{h^h)F{H-U)\ ' 

where the function F{t) is defined by 

F{t) = 1 - e''"(*+*^) . 

In the thermodynamic limit i — > 00, we expand Eq. ( |3^ ) 
in powers of eo- Then Eq. ( |l5|) can be brought into the 
following form 



n{uj) = -2n2e~'^/Scos2(7r/25) / dt 



00 /"OC 



(-^-^a)l+V.7o J, 
(Tir2)i/f 



[(— Ti — r2 — t — ia){—Ti — t — ia){—T2 — t — za)]^/ 9 ' 

I 



where we used the fact that the LL field operators sat- 
isfy the braiding relations, (|l9|), with the braiding phase 
function 



ip{t) = - sign(i) . 
9 



(38) 



Now we can perform the r-integrations using the integral 
representation 



r(i/g) Jo 



X e 



— i(ri+T2— t)p 



Then 



.(w) = 2i-f^e-"'3 cos^{iT/2g) 



(39) 



X / dt 



{-t-iaY+^/s 



X r2(l/5 + l) / dpp^'^^^eP*^^{llg,0]-i{p + iv)t) . 
Jo 



where F and ^ are the gamma function and the Tri- 
comi confluent hyper-geometric function, respectively 
psf . Deforming the ^-integration contour from the real 
axis to the contour around the branch cut of the function 
^ we arrive at a more convenient representation for the 
TED: 

rV-uj pi/g-i 
n^{oj)^A{g) dE Fv{E + u), (40) 

Jo (E + Luy 

where the spectral function is given by 

/•OC 

Fv{p) = 2lme-'''^s / ^-i/g-ip-{v-p)i 

J a 

x^\l/g,0,-p^ + ia), (41) 

and the pre-factor is 

A(,)=F^(l/,+ l)7V/^SS^ 



27ralT{l/g) ' 

In the limit of small energies, uj/V <C 1, the function ^ 
is nearly constant and the limiting form of the spectral 
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function can easily be established: 

27r 



Fv(p) 



r3(l + l/.g) 



Therefore we obtain the following asymptotic form for 
the TED in vicinity of the Fermi edge 



(42) 



where A = 1 / g—2 and C2 is a numerical pre-factor regular 
at A = -1. 

We investigated the behaviour of the TED near the up- 
per threshold uj = V hy numerical evaluation of Eq. (40). 
The result is given by another power law: 



n{uj)^ {V-ujY , v = X + 2. 



(43) 



In the limit of weak interactions, 5 — > 1, the exponent 
A approaches the perturbative result. Eq. ( p^ ) can be 
regarded as consisting of two factors. The first factor 
is the universal l/uj divergence inherent to all interact- 
ing ID systems. The second factor reflects power-law 
renormalisations occurring in the LLs and, in particu- 
lar, contains the local density of states (LDOS) of the 
primary electron. The latter object is suppressed for re- 
pulsive interactions [g < 1). So also the TED singularity 
is suppressed. At gc = 1/2 the LDOS suppression ef- 
fectively wins over and the singularity disappears. It is 
noteworthy that in the case when the interaction con- 
stant is smaller than this critical value, the TED has a 
maximum between the upper and lower thresholds uj — 
and cj = y, as it is vanishing towards both limits. 

In real systems the tunnelling amplitude 7 is small. 
The secondary current is already calculated in the next- 
to-leading (compared to to the primary current) approx- 
imation in 7. Still, it is a valid question to ask what 
happens at higher orders in the tunneling amplitude. 
For the case when the TED is divergent at the Fermi 
level (A < 0), a small energy scale lu* may emerge be- 
low which these higher order contributions become im- 
portant. We do not currently have full answer to this 
question. Our preliminary calculations indicate, how- 
ever, that the most divergent 4-th order diagram is the 
one for which Fig. ^ acts as self-energy. This would lead 
to the estimate ui* ~ ^7"^/-^. This issue deserves further 
investigation. 

Thus far we have discussed the spinless LL. The above 
results can be straightforwardly generalised to spinful 
systems such as quantum wires. Due to the spin-charge 
separation the LDOS at the end of the wire is then given 
by the same expression as in the spinless case with 1/g 
substituted by (l/ffc + ^/gs)/^, where gc,s are the inter- 
action constants in the charge and spin sectors, respec- 
tively p^ . Since the tunnelling processes conserve spin, 
the LDOS exponent is simply carried over to the TED, 
so that A in Eq. (E2|) should be changed to: 



(44) 



The relation between the lower and upper threshold ex- 
ponents, Eq. (^), is stiU vahd. When the spin SU(2) 
symmetry is preserved, it forces gs = 1. The residual 
spin-backscattering interaction renormalises to zero un- 
der renormalisation group transformations (see for 
a recent review on ID physics and further references). 
However, the presence of an irrelevant operator should 
probably cause multiplicative logarithmic corrections to 
our power-law formula (^). We have not attempted to 
calculate these correction though we expect that this can 
be done by standard methods p6|. 



VI. TUNNELLING FROM 2D SYSTEMS: PURE 
LIMIT 

In view of applications to MWNTs, we now consider 
tunnelling from a 2D electron system in tip geometry. 
Theoretically, this is a more complicated situation than 
in ID. In some cases, relevant non-perturbative tech- 
niques exist, in other cases we shall present perturbative 
results. In this Section we discuss interacting 2D electron 
systems free from impurities (pure limit). 

From the mathematical point of view we still have to 
calculate the same diagram [(a) in Fig. To accom- 
plish this task in 2D it is convenient to change over to 
the momentum representation. The host Hamiltonian is 
H[^]^Ho[^] + H,[iP], where 



Horn = / 



2tt- 



(45) 



describes free 2D electron gas with dispersion relation 
^{pj and V'p is the Fourier transform of the electron field 
operator. The interaction term has the standard form. 
In real space: 



dfdr' i;\f)ip\r')Uoi\f - r'|)V'(r')V'(^) , (46) 



where Uo{r) is the interaction potential to be specified 
later. The relevant diagram is given by the following 
expression: 



d'^q'i f d^q2 f dVl f d^p 



A ^ A = (1/g, + l/g,)/2 



(27r)2 J (27r)2 J 2t: J (27r)2 

X gQ~{^,P- oDOa'^i^ -^^P) (47) 
X Uo{qi)Tl{n; ql, q'2)Uo{q2)gi^+ {uj,p - g^) 

Here Uo{q) is the Fourier transform of the interaction 
potential, H(il; (71, ^72) stands for the inner polarisation 
bubble as in Fig. ^ (b) with vertex indices as in (a): 

^ ^ f de f d^p 

H(r!; 91,92) ^ j ^ j j2^9o {f ,P + <li) 

X G^+{e)g++{e,p + q2)g^+{e + n,p). 

Without the tunnelling insertions, this would be a 
Keldysh analogy to the standard 2D Lindhard functions 
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p7[ . We define p = Pujfi^ with a unit 2D vector ri0 in 
direction parametrised by the angle (j). The quantity 
gives the on-shell magnitude of the momentum p and, 
for a quadratic dispersion relation, satisfies the following 
equation, 



p2 

2m 



2m 



(48) 



where is the Fermi momentum and m is the effective 



We change over from the momentum integration to the 
energy-angle integration in the standard manner: 



(27r)2 



d<t>, 



where p(^) is the density of states in the host. The po- 
larisation bubble can then be written in the form 



-V Jo 



r 



where pc stands for the constant density of states in the 
lead. The diagonal bare Keldysh Green's functions are 



{ + + )/■ D 

9o '(e,Pe+nn^ 



-p. 



2m 



In this expression we have taken into account the fact 
that e < if we are only interested in oj > part of g ^ . 
The off-diagonal Keldysh function is 

go'^iuj - ri,p) = i2TinF{p) 5{uj - - ^p) , 

where np{p) is the Fermi momentum distribution func- 
tion. Using these expressions and relation (H) between 
the Green's function and the TED we obtain a formula 
similar to Eq. (^J) for the ID case, namely 

n>(w) = — / dnp{uj-n) [ dep{e + n) (49) 



where the function F is defined by 

F,(17,0) = — dqqUoiq) d0 

X [n + —P,+nqcos{e^(f,/2) + ^q^+iO]-'^ 
m Zm 

X [n + —P,+nq cos{e + (t>/2) - + iQ]-^ . 

TO 2m 

To proceed, we now need to specify the interaction 
potential. It is natural to start with a contact inter- 
action term, U{r — r') = U{)5{r — r') or Uo{q) — Ua 
in the momentum space. However, an unexpected tech- 
nical difficulty arises. If we follow the standard prac- 
tice and linearise the electron spectrum in the vicinity 



of the Fermi surface, then the momentum integration in 
Eq. (|5^) can be easily done but the remaining angle in- 
tegration diverges at ip = ±7r. This latter singularity 
comes from large momenta in the previous momentum 
integral and corresponds to the special case of hot hole 
back-scattering. To obtain a finite result for the TED 
one either needs to take into account the decay of Uo{q) 
for large q or the non-linearity of the electron spectrum. 
Accordingly, we present two different calculations in the 
rest of this Section. 

(i) Contact potential: Uo{q) — Uq. 

First we scale the integration variable q —f flq and 
introduce dimensionless quantities 71^2 = Pe±Q/PF and 
K = n/AEp. Then, 



Un 



00 /'27T 

dqq / d9 
Jo 
2 



X [l + 7iqcos(6'-(/«/2) + K(7^+iO]"^ 
X [l + ^2qcos{e + (f>/2)- nq"^ +iO]-K 

Notice that 7's approach unity whereas k linearly tends 
to zero in the limit of small fl. Furthermore, if we set n = 
0, the angle integration is still divergent. Therefore we 
set 71,2 = 1 and keep k finite. As a next step we expand 
the integrand in Eq. dSO) using the standard formula [B8| , 



1 



= V- 



1 



f{x)+iO fix) 



- m5{f{x)) , 



(50) 



where V denotes the principal value. This expansion 
produces three contributions: the first one contains the 
product of two principal parts and remains regular for 
57 ^ 0, the second one contains products of one prin- 
cipal part and one delta function and is identically zero 
in the low energy regime, and the third contribution, ^ 
containing a product of two delta functions is responsible 
for the divergency or the angle integral. Setting vp = ^ 
for the rest of this Section, we write 



13 



= — ^/ dqq dd6[l + qcos{9-(l)/2) + Kq^]S[l + qcos{e + (P/2)-Kq^] 



lo Jo 

The evaluation of the angle integration leaves us with 



47r2 



dqq[q^ - {I + nq^f] ^^'^6[l- (1 + nq^) coscf) + aincpy/ q^ - (1 + f^q^Y - 



r 



The argument of the delta function is zero for 
isin^(0/2)( ' 



ql 



1 - Jl - 16- 



Therefore the remaining momentum integration yields 



27r2 



smffl + cos ffl sm 



sir 



sin (f) — 16k' 



4k sin c 



Plugging this expression into Eq. (|4^) and taking into 
account that for small (and hence small n) the main 
contribution to the angle integration originates for (f> close 
to ±7r one finds for the following limiting form for the 
TED: 

where = p(0). 

(ii) Yukawa potential: Uo{q) = 4776^/ q^ + A^. 

Here 1/A is the screening length. As U[){q) now tends 
to zero for large q (which is, of course, always the case 
in real systems), we can linearise the dispersion relation. 
Then the function F{Vl, 0) becomes independent of the 
energy variable e and simplifies considerably: 



1 /■°° 
F{n,^) = —J^ dqqUo 



mi{q,4>), 



(51) 



where q has again been scaled by Q and the function 
I{q, 4>) is defined by 



d0 



(1 + gcos(6'- 0/2) + iO) 



X (l + qcos(6' + (?!)/2)+iO) 



This integral can be calculated and is real 



[1 cos2(</>/2)]yr^ 

for q < I but contains an imaginary part 

/(g,0) = ■ 



[l-g2 cos2(,/,/2)]V?^ 



(52) 



(53) 



(54) 



y^23icos(0/2) 



,5[9-l/cos(0/2)]. 



for q > 1, as incoherent particle production takes place in 
the latter regime. Using Eqs. ( |5^) and (^J) we compute 
the momentum integral in Eq. (|5l|) which yields 



F{n, 



2e^ 



X arccot 



1(0/2)1^172 + A2 cos2(0/2) 
^n2_^;^2 cos2(0/2) 

A|sin(0/2)| 



(55) 



Unfortunately we were not able to perform the last re- 
maining angle integration in Eq. ( |49| ) in a closed form. 
However, to analyse the TED close to the Fermi surface 
— > one merely requires the knowledge of \F(n, 0)|2 
for small fl^s. The latter can be easily read off Eq. ( [Sq ) 
for f2 <C A. There are two different regions in the param- 
eter space (0,11): for most angles, |0 ± 7r| > 211/A, we 
find 



A2 



sm 



whereas for large scattering angles |0±7r| < 2fl/X one 
obtains a singular in O behaviour, 

|F(f7,0)p 



f72 



The latter region, physically corresponding to hot hole 
back-scattering processes, is responsible for the leading 
contribution to the remaining integrals. With logarith- 
mic accuracy we thus obtain the following threshold be- 
haviour of the TED, 



,27re^ 



■In 



max(F, A) 



(56) 



Clearly the same logarithmic singularity will persist for 
other functional forms of Uo{q) as long as it decays for 
large q. (We note that in the purely Coulomb case, A = 0, 
an additional divergence at small q is introduced. Cal- 
culation shows that all the divergences cancel exactly in 
this special case with the second order result for the TED 
being identically zero for all a; > 0. The physical meaning 
of this observation escapes us.) 

The net result is that we find a logarithmic singular- 
ity in the TED for 2D electron systems at the second 
order in the interaction. Such behaviour is in apparent 
contradiction with Ref. Q where a stronger singularity 
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was found. These differences may be due to the fact that 
Gadzuk and Plummer analysed a different set-up with 
a three-dimensional host and a 2D emitting surface (as 
opposed to our set-up with a 2D host and tip emission), 
or their stronger singularity may be an artifact of the 
low-density approximation. This issue as well as the role 
of higher-order scattering processes deserve more inten- 
sive study and will be discussed elsewhere. In view of 
applications to MWNTs, however, there are two more 
immediate questions: the role of the dynamical Coulomb 
screening and of the disorder potential. Both require a 
non-perturbative approach and will be discussed in the 
next Section. 



VII. TUNNELLING FROM DISORDERED 2D 
SYSTEMS 



Contrary to the SWNTs, MWNTs can be described 
by the LL theory only in some special cases when the 
radii of outer shells are not too large Otherwise 
their physics is best described in terms of a disordered 2D 
electron liquid with Coulomb interaction. Such systems 
have been extensively studied in the 80 's using non- linear 
sigma model and renormalisation group ^ . Re- 

cently they attracted more attention due to the dis- 
covery of a possible metal-insulator transition in MOS- 
FETs , which stimulated further theoretical research. 
Kamenev and Andreev (KA) |Q have recently exam- 
ined the Keldysh non-linear sigma model for the 
Coulomb case (long-range interaction). In this Section 
we apply their method to obtain the threshold behaviour 
of the TED for tunnelling from a 2D interacting disor- 
dered metal. 

The Hamiltonian of the previous Section, composed 
from (45) and ( |4^ ) with Uo{q) — Ane'^/q, should now be 
supplemented by the disorder term 



Hd[xP] = / drV{r}ip^r}ip{r} 



(57) 



with a (5-correlated Gaussian disorder potential V{f). 

The formal procedure of deriving the Keldysh sigma 
model is described in the relevant literature (Q, see 
also Ql) in detail. Therefore we only give a brief outline 
here. One starts with the Keldysh S'-matrix. (This ob- 
ject is identically unity but it becomes a generating func- 
tional upon introducing auxiliary fields. If one wants to 
start with a diagonal bare Keldysh function, the contour 
is different from that in Fig. ^, see explanation in Refs. 
1^ , Q). One can then integrate the disorder out ob- 
taining a non-local in time four-fermion interaction term. 
The action is made quadratic in fermions by applying the 
Hubbard- Stratonovich transformation to the latter term 
as well as to the Coulomb interaction term. Hence there 
are two decoupling fields, Q and $ in the notation of 
Ref. iQ. The fermions can now be formally integrated 
out resulting in the following expression for the generat- 



ing functional Q]: 

{Z) = J I?$e'Tr(*-c/-Vi*) J -pQ^^siQM (gg) 

with the sigma-model action (after Keldysh rotation) 



Tr ln[Go"^ 



2n^' 



+ <jyc.r + C] 



Here tq is the elastic mean free time and [/q stands for the 
unscreened interaction potential. The decoupling field Q 
depends on two time variables and is a 2x2 matrix in 
the Keldysh space. The decoupling field $ is a Keldysh 
doublet {(t>i,(t)2)^ ■ The vertex matrices are defined as 
follows: 



1 
1 



7 = (Ti = 



1 

1 



Doing the functional variation of the generating func- 
tional, (^8[), with respect to the auxiliary field C and set- 
ting the latter to zero one obtains the complete set of 
single-particle Green's functions: 



C=o 



' 



where the function = (G^ - G^)(l - 2n{E)) is re- 
lated to the single-particle energy distribution function 
n{E) and G^^"^-' are the retarded and advanced Green's 
functions, respectively. 

We need to calculate of the four-point correlation func- 
tion /C^^ (t2) ^3 5^1) ^4)- (A calculation of a different four- 
point function for a non-equilibrium noise problem has 
recently been done within a similar framework in pTf.) 
This can be achieved by double variation of the generat- 
ing functional with respect to the auxiliary field C,: 

IC+l{h,h;h,U)= J I?$e*Tr(*-v^-i*) J I?ge^S[Q,$] 

X [Wih,U) ® Wih, t2) + Wih,t2) ® Wits, UM59) 

where 



W{ti,t2) 



A functional integral of this type with a single W func- 
tion yields the Green's function matrix Q and was dis- 
cussed in [Q. A product of two and more W functions 
produces a complicated object, which is a tensor prod- 
uct containing all possible time orderings. The operation 
of extracting the particular component corresponding to 
/C^^ we denote by [...],. At this stage its exact defini- 
tion is unimportant and we postpone it to the end of the 
Section. 
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We now evaluate the functional integral with the prod- 
uct of two PF-functions in the saddle point approximation 
following KA's approach Q. In this case 



W{t, t') = -iTiv e''^"^*)^" K{t - t') e'fc°(*')'r° 



(60) 



where K{t) denotes the mean-field non-interacting 
Green's function matrix in equilibrium (corresponding to 
non-crossing disorder diagrams ) and it) is a linear 
functional of $. The corresponding correlation matrix 
with respect to averaging over the (/)-fields is given by 



{ka{q,u))ki3{-q,-uj))^ = -Va/3(g,w), 



(61) 



V«(<z,w) V«(g,cc>) 
V^(q,w) 



v«(^)(<z,^) 



1 

[dq^ ^ iivY \_Vo dq'^^iu) 
= ns(a;)(V«(g,a;)-V^((Z,^)), 



i^dq'^ 



(62) 



where nB{<^) is the Bose distribution function and d is the 
electron diffusion constant. Using the saddle point rep- 
resentation in Eq. (^0|) the product of two VF-functions 
can be re-written as 

{[W{hM)®W{t^M)]:)^ 
1 ^ , 



4 ^ 

Mi 



(63) 



The operation of taking the correct time ordering is now 
separated from of the averaging over the field <&, which 
only affects the objects 



„i(fei(t) + fc2(t)) 



„l(fcl(t)-fc2(t)) 



= s{t) + ^is{t) = e"('=iW+'3'=^(*» + ^e^^C^iW-Z^fe^W) . 

The product of four p^-operators is a sum of 16 products 
of different s(ii)-operators. The latter objects can be 
evaluated using the correlation matrix, Eq. (|6l|), 

/g4ai(fel(tl)+/5lfe2(tl))g-M4(fcl(t4)+/54fe2(t4)) 



X e 



ia3(fcl(t3)-|-/33fc2(t3))„-M2(fel(t2)+/32fe2(t2)) 



exp 



{ 



aia3Vi3(/3i,/33) - aiQ;4Vi4(/3i, /34) 

- Q!4a3V43(/34,/33) + a4a2V42(/34,/?2) (64) 

- a3a2V32(/33,/32) - aia2Vi2{Pu P2) + 2Vo(/3i, } , 
where we use the following definition 

Vy(A, /3,) = V^(<, -tj) + p.V'^iU - tj) 

+/3,V^(t,-i,). (65) 
Taken separately these functions are divergent. Never- 
theless Eq. (|6^), being rewritten as a function of differ- 
ences Vij{(3i,(3j) — Vii{(3i,(ii), can be made convergent. 
In this case a special selection rule has to be fulfilled, 



aia^yoiPiiPs.) - aiQ;4Vo(/3i,/34) 
a4a3Vo(/34,/33) + a4Q!2Vo(/34, /?2) 
a3a2Vo(/33,/32) - aiQ!2Vo(/3i, /32) 



0, 



(66) 



otherwise (|6J) is divergent. Using Eqs. (|6^) one can cal- 
culate the asymptotic behaviour of Vy (/?i, (3j)—Vii{(3i, Pi) 
for large time differences (to calculate the limiting Fermi 
edge asymptotics of the TED this knowledge is sufficient 
as we shall see shortly). Fortunately, the limiting form 
turns out to be independent of the arrangement of the 
/3-indices, with the result 



Vy(A,/3,)-V,,(/3„/3,) 



i In^ 

47r/ 



f\U-t,) 



is then 



where / = 2TTe^vd. The last average in Eq. (^ 
given by 

(Pmi (^1 ) (Pm4 {U ) ) *Pm3 (^3 ) (Pm4 2 ) ) * ) $ 
[](l+/ij) I (S(ti)s*(i4)s(t3)s*(t2))*. 



Obviously, the above is nonzero only if for all i — 
+1. That means that we can set all a/s equal as well. 
Such choice automatically fulfils the selection rule, (|6^). 
Therefore the full four-point correlation function is given 
by 



lCj^t{ti,t2\h,U) w F(ti,i2;i3,i4)exp y - -^-j ^ 



In" 



I\h-u){ti-h){h-t2){h-t2) 



d'^{ti-h){U-t2) 



}) 



(67) 



r 



with 



Y{hM-M, U) = [{-i^k{hM)i^) ® ii^Mh, t2)j')] . 



The latter expression is basically the non-interacting 
four-point correlation function with special time order- 
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ing and is equal to G"' {ti — t4)G"' (^3 — t2)- The sec- 
ond term in Eq. ( ^ ) can be calculated in a similar man- 
ner resulting in the same, up to an interchange of time 
variables, exponential factor as in Eq. ( |67| ) and an ap- 
propriate pre- factor. The pre- factors can in principle be 
re-exponentiated but this would only results in logarith- 
mic corrections which can be neglected in comparison to 
the main disorder part. In order to use formula (|l5|), we 
substitute ti = ti, ^3 = 0, t2 —t + Ti and t4 = < + ti +T2. 
We are interested in the low-energy regime u)/V ^ 1. 
Therefore the lower limit for the time integrations can 
be set to 1/D (D being the conductance band- width) if 
it causes divergences and to zero otherwise. Thus we 
obtain 



n>(w) 



dri 
^0 



X exp 



8^/ 



In" 



/*(ri+r2) 



Estimating the r-integrals for small cu yields the final 
result for the TED in vicinity of the Fermi edge: 



n^{uj) exp 



In" 



(68) 



Thus we find that the combined effect of the disorder and 
the Coulomb interaction is to suppress the TED towards 
the Fermi surface. This result can be understood in terms 
of the well known Altshuler-Aronov-Lee anomaly [|8| in 
the density of states of the primary electrons. Indeed, one 
of the important results of the KA's work was to show 
that the negative logarithmic correction to the density 
of states can be exponentiated (and becomes a double- 
log due to the long-range Coulomb forces). Our formula 
(pSi ) achieves a similar exponentiation for the TED. As in 
pjl, this is only valid above certain energy scale e* , due 
to the on-set of the fluctuation effects (corrections to the 
saddle-point approximation) . In MWNTs the situation is 
further complicated by the dimensional cross-over effect; 
the corresponding energy scale e* ~ d/L^ [L being the 
tube's circumference) was estimated in see also |^ 
for further discussion. 



VIII. FIELD EMISSION FROM CARBON 
NANOTUBES 

Field emission setup is relatively simple. An electric 
field of strength F is applied between a metallic tip of the 
emitter and a counter-electrode which are usually placed 
in a vacuum. Superposition of the confinement potential, 
the image potential and the electrostatic potential due to 
the external electric field leads to a formation of a barrier 
with a finite width allowing the electrons to tunnel out of 
the emitter even at zero temperatures, Fig. ^. Therefore 
the FE process can be regarded as tunnelling into vacuum 
through a triangularly-shaped barrier whose upper part 
is rounded off by the image potential. However, at low 




F=0 



FIG. 6: Schematic representation of the field emission setup. 
Here F is the applied electric field, W is the work function 
and D is the conduction band- width. The dashed line repre- 
sents the confining potential without field and the dotted line 
stands for the image potential. 



temperatures the larger part of tunnelling electrons will 
have energies close to the Fermi energy, so the exact form 
of the barrier is unimportant. 

Quantities of interest are the energy-resolved current 
j{Lo) and the total current J. If the tunnelling amplitude 
is small (which it is for all real systems), then j(a;) is 
proportional to the probability n(w) for an electron to 
have energy to the energy dependent transmission co- 
efficient I?(w), and to a factor F responsible for the tip 
geometry: 



j{uj) = FV{uj)n{Lo) . 



(69) 



In the noninteracting case and neglecting higher-order 
tunnelling processes, n{uj) coincides with the Fermi dis- 
tribution function multiplied by the LDOS. The (quasi- 
classical) transmission probability for a triangular barrier 
is given by Q : 



ufl'^/ifiF) 



(70) 



where m is the electron mass and W is the work function. 
As discussed e. g. in Ref. Q, the transmission coefficient 
is only slightly affected by the image potential and for 
the field strength F <^ {W — EpY/e^ the barrier can 
be regarded as strictly triangular. At zero temperatures, 
the emerging (primary) spectrum has a sharp threshold 
at the Fermi energy (no particles above Ep) and is es- 
sentially constant in its vicinity. 

The general picture does not change much even in the 
case of interacting electrons. For the LLs n(w) is propor- 
tional to the relevant LDOS, which is known ^8, Q to 
be 

n{Lo) = ei-Lo)\Lu\^^s-i /a^jji/9Y{i/g) . (71) 



Plugging this into Eq. (|69|) one observes that at this low- 
est order in the tunnelling the TED above the Fermi en- 
ergy is still zero, even for the strongly correlated LL. 
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Below the Fermi energy the TED has a power-law singu- 
larity. Expanding the transmission coefficient in powers 
of uj/W (work function being the largest energy scale), 
and integrating over all energies we establish the Fowler- 
Nordheim (FN) formula for LLs: 



J = 



1/2S 



4fc 



1/2 



exp 



3F 



-M/3/2 ^ (72) 



relating the full current to the electric field's strength 

The generalisation of these results for SWNTs is rather 
straightforward. Indeed, these systems are known to 
be described as four-channel LLs. Three channels (f)c~ 
(charge-flavour), (j)s+ (total spin), cjjs^ (spin-flavour) are 
non-interacting. The fourth channel (f)c+ (total charge, 
or the plasmon mode) possesses the LL parameter K = 
(1 -|-4[/o/7r)~^/^, where Uq is the zero Fourier component 
of the screened Coulomb potential. Note that though we 
now have four channels the field-operator actually fac- 
torises as |l|, |l|, |l|] 



iP ~ exp{i0c+/(2\/X) + ii(t)c- + (t)s+ + 0s-)/2} . (73) 

Just as in the case of the spinless LL, the interaction 
constant K is included into the rescaled Bose fields and 
disappears from the Hamiltonian, which is given by 



H 



1 

47r 



E 

(5 — C,6*J — ± 



dx {dx(j>Sj] 



(74) 



All the correlation functions also factorise. Hence the 
results for the LDOS [2^ , |3^ as well as the results of 
Section 0, Eqs. (||) and are still valid for SWNTs 
given that the substitution 



(isT-i +3)/4. 



(75) 



is made. 

The typical value for K in SWNTs lies near 0.2 
(0.15 to 0.3), which fixes the effective q to « 0.5 
||, |l|, |l|, |l6). As a result, Eq. (|7|) is reason- 
ably well approximated by the classical FN law. This 
fact offers one possible explanation as to why the exper- 
imental data of Refs. pl| can be relatively well fitted 
by the conventional FN curve. In addition, all the ex- 
periments currently available were made on SWNT films 
where the nanotubes build a network with essentially a 
random 2D geometry. This case, when extra complica- 
tions are bound to arise, is beyond the scope of this pa- 
per. To straightforwardly reveal the correlation effects in 
the primary current, measurement on a single SWNT are 
more appropriate. To our knowledge, this has not been 
done yet. 

Let us now discuss the secondary current. The analy- 
sis of the secondary effects can also be made in terms of 
tunnelling into vacuum. Thereby we should bear in mind 




that the theory presented in Section [II assumes that 



FIG. 7; A sketch of the energy resolved current in a field 
emission process from LLs. An additional Lorentz peak due 
to tunnelling through a localised level (if exists) is depicted 
in the inset. 



the tunnelling amplitude is energy independent. There- 
fore, since the effective applied voltage is large in the EE 
setup (one has to send the chemical potential on the right 
electrode to minus infinity), some integrals inevitably di- 
verge. To overcome this difficulty and to made the theory 
more quantitative, we now recall that there are two can- 
didates that can be used as the effective voltage. The 
first one is the characteristic energy scale in the expo- 
nential transmission coefficient A = F/2{kpWy^'^ , see 
Eq. (|7^). The second one is the width of the emitter 
conductance band D. Therefore the high-energy cutoff 
playing the role of V in the FE setup is either D or A, 
whichever is smaller. We shall use D through the rest 
of the paper. With this modification all results of the 
previous Sections for the TED n{uj) are valid with the 
understanding that the tunnelling probability is propor- 
tional to the transmission coefficient of the barrier at the 
Fermi energy, 7^ ~ 'D{Ep). According to Eqs. (Q) and 
(^) the secondary current is proportional to the TED. 
For SWNT the TED possesses all the singularities char- 
acteristic for LLs, which we have already discussed. One 
thing worth noting is that the critical value of the cou- 
pling now is Kc = 1/5, which is actually within the ex- 
perimental range. Therefore we do not make a specific 
prediction regarding the character of the singularity (di- 
vergent versus convergent TED at the Fermi edge). In- 
stead we think that both the divergent and the vanishing 
TEDs can be observed depending on the experimental 
setup. Unfortunately, in most recent experiments, both 
on SWNTs and MWNTs, only the total current is mea- 
sured. Where the energy-resolved measurements were 
actually made ||2^, |5^, |5^, in all cases but one, the tem- 
peratures were too high for the secondary effects to be 
visible. To our knowledge, so far only the experiments 
of Rcf. contain high-energy tails which can be at- 
tributed to secondary emission. However, the quality of 
the presented data hinders us from attempting to actu- 
ally fit the curves. 
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FIG. 8: A schematic representation of the electron tunnelling 
through a localised state d. 

We conclude this Section by summarising our scenario 
for the energy resolved current in the FE from LLs (and 
SWNTs). Well below the Fermi energy j{uj) is governed 
by the exponential growth of the transmission coefficient, 
j{uj) exp{uj/A), see Fig. ^ as is the case for any metal- 
lic emitter. Nearing the Fermi level, at energies given by 
iJp « (1 — l/g)A, the non-linear LL LDOS effects win 
over the exponential growth resulting in a maximum in 
the current profile Further towards the Fermi en- 

ergy it decreases according to the power-law: ~ 
The edge behaviour right above Ep is given by ~ uj^^^'^. 
This changes shape depending on whether g is smaller or 
larger than the critical coupling gc- In the case of weak 
interactions, 1/2 < g < 1, there is a singularity, while for 
strong correlations, < g < 1/2 this singular behaviour 
is suppressed and there is a power-law approach to zero. 
In the latter case j{uj) acquires an additional maximum 
at a cross-over energy E* . For small applied fields, there 
is an upper threshold at energy imii(D,A) where the en- 
ergy resolved current behaves as ~ {mm{D,A) — ivY^^ . 

IX. TUNNELLING VIA LOCALISED STATES 

Recent luminescence spectra measurements in FE ex- 
periments on carbon nanotubes ||2^, |2^ suggest that, at 
least in some cases, the emission process cannot be under- 
stood simply in terms of tunnelling into vacuum. They 
rather seem to be compatible with a model where the 
emitted electron tunnels out of the system through one 
or more localised states at the tip of the emitter [g^ ^ . 
In this Section, we generalise our formalism in order to 
analyse such a multiple stage tunnelling. 

We start with the Hamiltonian of the system. Contrary 
to Eq. (|l|), the transfer of an electron between the host 
and the lead occurs in two stages. At the first stage, the 
electron populates a localised level with energy —A, see 
Fig. ||. Let the tunnelling amplitude for this process be 
7o. At the second stage the electron tunnels from the 
localised state into the lead, 

H2 = + Ho[c] + 70 [^-(0)^^ + dtV(O)] 

+ 7 [d^c{0) + c\0)d] - AdU , (76) 



where we retain the notation 7 for the tunnelling ampli- 
tude of the second process and d'^ , d are the creation and 
annihilation operators of the localised state, respectively. 
We assume that < —A < —V. This is not restrictive 
but reasonable. Indeed, the first inequality reflects the 
fact that the electric field applied to the emitter causes 
an effective lowering of energy of the states close to the 
tip. The second inequality is a trivial one because if the 
energy of the state d lied below both Fermi levels it would 
have always been populated. For simplicity, we neglect 
any electrostatic interaction between the localised state 
and the leads and the Hubbard term for d. Like in the 
case of the simple tunnelling, the lead is assumed to be 
uncorrelated. 

The problem of a localised state hybridised with a non- 
interacting continuum is, of course, exactly solvable [Q. 
Therefore one way to proceed would be to eliminate the 7 
hopping term from the Hamiltonian exactly and then ap- 
ply perturbation theory in 79. However, in the FE setup 
A is expected to be much smaller than the conductance 
band width D. This implies that 70 ^ 7. Therefore 
we take an alternative route and first calculate the TED 
of the electrons on the level d as a function of g '"(w) 
in the host ignoring tunnelling to the lead. (We did the 
alternative calculation as well and, at a given order of 
perturbation theory, obtained exactly the same results.) 
To proceed we define the Keldysh Green's function of the 
localised level: 

D{t) = -i{Tc[d{t)dH0)Sc]) , (77) 

where the 5-matrix includes only the tunnelling between 
the host and the localised level, 

Sc = Tc exp (^-ijo dt [4>^ (0, t)d{t) + S (i)^(O, t)] ^ . 

Expanding the Green's function ( |77| ) and function 
{Tc[d{t)il}'^ {Q^t')Sc\) in powers of 70 one obtains a set of 
coupled equations for them. Eliminating the latter ob- 
ject one then obtains the following Dyson-type equation 
for the Keldysh Green's function: 

D[t, t') = Do{t, t') + [ dt" D{t, t")K{t", t') , (78) 
Jc 

where the kernel is given by 

^{t'\t') = 72 / dtg^{t",t)D^{t,t') . (79) 
Jc 

Here go{t,t') and Do{t,t') denote the Green's functions 
of the host and the localised level in the absence the hop- 
ping term. Disentangling the Keldysh indices and Fourier 
transforming we obtain the following set of equations: 

D-+{lu) = D^+{uj) + D—{u)k-+{uj) - D-+{lu)k++{uj) , 
D—{uj) = Dq-{uj) + D—{uj)k--{uj) - D-+{uj)k+-{uj) . 

The solution is given by 
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D-+ico) 



I 



(80) 



This equation can be further simpUfied if one takes into 
account that Jie[gQ~ {u})]/\Dq~ {lo)]"^ vanishes for our sys- 
tem. Then the TED on the locahsed level is proportional 
to 



(C. + A)2+74|^-(^)|2 



(81) 



In the last relation we omitted the subscript of g{uj) 
functions. Strictly speaking this constitutes an approx- 
imation. The reason is that since the tunnelling onto 
the localised level affects the Green's function of the host 
they have to be calculated self-consistently. However, 
the renormalisations occurring close to the Fermi energy 
should not strongly affect the shape of the localised state 
and vice versa. That is unless the localised level becomes 
resonant (A — s- 0), which case we shall not consider in 
this paper. (We note though that the resonant level case 
can be approached via a mapping onto the Kondo prob- 
lem |5^.) For a non-interacting host \g^^{Lu)\'^ is simply 
a constant equal to 1/ivp. 

the TED in the non- 



As we know from Section III 



interacting lead, is proportional to D '^(lu) and the tun- 
nelling probability 7^: 



n{uj) 



(82) 



We now briefly discuss the picture for a non-interacting 
emitter. According to Eq. ( pl| ) there is a Lorentzian 
peak a,t Lu — — A with height 4v^g *"(— A)/7o and width 
Jq/2vp. At the second order in 7^, there still is a sharp 
threshold at the Fermi energy. These results for non- 
interacting emitters are, of course, known (see e. g. [^). 
Here we have merely re-derived them using the Keldysh 
formalism. 

Let us now turn to the open question of how the pres- 
ence of a localised state would modify the FE process 
from a LL. For the primary current we can still use for- 
mula (|69|), where we have to substitute the TED of the 
electrons in the host by that of the electrons on the lo- 
calised state. The energy resolved current is then given 
by the following expression (we again set vp = 1), 



Jiuj) = -iTV{uj)D-+{uj) 



11/9-1 



X exp —Ak 



3F 



(83) 



In comparison to the situation without the localised level, 
there is an additional feature: the Lorentzian resonance 
at w = — A. To compute the total emitted current we 



re- write the energy integral in terms of the dimcnsionless 
variable ^ = lu/A and expand the exponent in powers 
of uj/W as the work function is large compared to other 
energy scales. The total current then is 



J 



a^D^l9Y\\lg\ 



exp 



3F 



1/9-1 



(^-AM)2 + (72/2^)2 



(84) 



This integral can be calculated exactly, resulting in the 
modified FN relation 



J 



A 



|w;i/9-ieU'r[l- l/g,iy]| 



applicable for arbitrary parameters. Here J is the total 
current in the absence of the localised state, see Eq. (7^), 
r stands for the (incomplete) gamma function and we 
have introduced a dimcnsionless quantity w, defined by 

^ = \ {illl'^ - A) . 

There two important limiting cases. 

(i) The electric field is weak, A/yl > 1, or A > 
F /2{kpWY^'^ . Then the resonance peak is far away from 
the Fermi edge in comparison with the characteristic de- 
cay scale of the transmission coefficient. The result of 
the integration in Eq. (^) is then essentially identical to 
the result for a system without a localised level Eq. ([7^), 
up to numerical pre-factors. Therefore we would not ex- 
pect that in this regime one can differentiate between the 
direct tunnelling and tunnelling via a localised state. 

(ii) The electric field is sufficiently strong, so that 
/S./A ^ 1. Now the distance from the resonance to 
the Fermi edge is much smaller than the characteristic 
energy scale A. If also 7^ <C A (which is to be ex- 
pected) the total current is still given by Eq. ([z^), apart 
from the overall pre-factor, but with a different expo- 
nent: 1/2(7 should now be substituted by l/2g— 1. This 
change can be important for the analysis of experimen- 
tal data on the SWNTs. For SWNTs in this regime, we 
expect the modified FN plot to be linear in the coordi- 
nates ln[J^/F(i/-^"5)/4] and l/F. For K w 0.2 the total 
current depends exponentially on the field strength (not 
a stretched exponential). 

The calculation of the secondary current is more deli- 
cate. In the case of the localised state the insertion in the 
diagram for g '"(w), see Figs. ^ and |[ should be mod- 
ified. It turns out that there is only one way to dress 
the Green's function Gq~ with the Green's functions of 
the localised level, G'(|~(w) ^ L»++(w)Gj-(w)L»— (w). 
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so that the diagram is non-zero above the Fermi edge. 
The time- and anti-time-ordered locahsed level Green's 
functions take the form 



D- 



-(++) 



1 



+ A ± z(ad + 7o/2) ■ 
The imaginary part in the denominators reflects the fact 
I 



that the natural width of the localised level ad is now 
widened by the hybridisation with the continuum of 
states in the host. Using the Fourier transform of the 
new insertion, which must now be substituted instead of 
the Fourier transform of Gq^{uj) (given in Eq. ( p^ by 
expi{oj — V)t/{t + ia)), we construct the localised level 
counterpart of Eq. (ITq) : 



n>(w) 



dte*("+^)*(e("''+^o)*r[i(arf + 7^ + i{A + D))] - e-("^+^o)*r[i(-ad - 7o + + D))]) 



2Ti{ad + %) J- 

pOQ nOO 

/ dri / rfT2e""("i+"^)/C(Ti,0;t + ri,< + ri 
Jo Jo 



r2) 



r 



(85) 



The asymptotic behaviour of the ny{uj) close the Fermi 
edge is not affected (apart from a numerical pre-factor) 
by the presence of the localised state. There is an im- 
portant difference - the upper threshold for the TED is 
now given by A instead of D. The reason is that in the 
presence of the localised state the upper energy limit of 
the hot hole is given by the energy of the d-level instead 
of the band- width D (or the applied voltage V). 

Thus, in the case of tunnelling via a localised state the 
TED obtains an additional feature, namely a Lorentzian 
peak below the Fermi edge, see the inset in Fig. |7[ Re- 
cent luminescence experiments of Ref. suggest a sim- 
ilar structure - two superimposed peaks in the spectrum 
of emitted light. The Lorentzian shape of the peaks was 
probably hidden by the intrinsic Gaussian broadening in- 
troduced by the apparatus (a Gaussian fit was adopted 
in Ref. [Q). The calculation of the emerging lumines- 
cence spectra is beyond the scope of this paper and will 
be discussed elsewhere. 



X. SUMMARY AND CONCLUSIONS 

In this paper, we studied the energy resolved current 
(or the TED) emitted from a tip of a correlated host ma- 
terial under out-of-equilibrium conditions. In the uncor- 
related case the TED has a sharp threshold at the Fermi 
edge and there are no particles with energies above it. In 
the presence of interactions there is a finite high-energy 
tail due to the hot hole relaxation process. For most 
systems, the TED is singular right above the Fermi edge. 

In the repulsive Luttinger liquid model, where the 
LDOS is renormalised by interactions and vanishes to- 
wards the Fermi edge, the TED below the Fermi edge 
shows a power-law behaviour with the LDOS exponent. 
Above the Fermi energy, we also find a power-law but 
with a new exponent. There is a critical coupling and 
the TED is divergent for weak interactions but conver- 
gent for strong interactions. In the latter case, the TED 



has a maximum above the Fermi level. Although the 
high-energy tail still exists in systems with local interac- 
tions confined to the tip, the TED is then a regular func- 
tion of the energy. Additional features in the TEDs arise 
in a situation when the tunnelling occurs in two stages 
via a localised state. In this case there is a Lorentzian 
peak centred at the energy of this localised level and of 
the width which depends on the hybridisation with the 
continuum states on both sides of the contact. All the 
singularities of the TED above the Fermi energy survive 
the introduction of the localised level. The singularities 
in high-energy tails are suppressed as soon as the disorder 
is introduced. We have calculated the TED of particles 
tunnelling from a correlated disordered 2D system. In 
this case the LDOS is exponentially suppressed and this 
suppression is carried over to the TED. Contrary, in the 
pure limit, we find a weak (logarithmic) singularity in 
the TED above the Fermi edge. 

We then specialised our results to the case of the 
EE from carbon nanotubes, both SWNTs and MWNTs. 
While they are in qualitative agreement with existing 
experiments, there is little data currently available to 
test our predictions for the secondary current in these 
systems. We hope that this work may encourage more 
measurements, especially at low temperatures (<C 300K). 
This could provide additional information about the na- 
ture of interactions in nanotubes. 

Throughout the paper we consistently worked at zero 
temperatures. All our results can be generalised to finite 
temperatures by standard methods (though this may not 
be of immediate interests until there is more experimental 
data). There are, of course, other open questions (like the 
role of higher orders in the tunnelling amplitude) and 
potentially interesting extensions of this work, which are 
detailed in the main text. 

Also interesting, though more difficult, could be exper- 
iments on tunnelling junctions if the TED in such setups 
were accessible e. g. by means of a tunnel probe Js^ . In 
this case both threshold asymptotics could reveal essen- 



21 



tial details about the correlations. Even if direct mea- 
surements of the TED are not possible, there is another 
possibility to identify the secondary effects. Since the hot 
hole decay is a relaxation process, an electron-hole recom- 
bination can occur with an irradiation of a photon. Such 
a current-driven luminescence effect could be observable 
in quantum dot systems, where the corresponding experi- 
mental techniques became widespread in the recent years 
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